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1 Introduction 



The problem of a quantum impurity in 3D electrons has been a fascinating topic in 
condensed matter physics. It is related to many nontrivial phenemena such as the 
Kondo problem. The low temperature behaviours of the Kondo and the Anderson 
models are rigorously obtained by the Bethe ansatz method. |jl], |^ Recently, some 
aspects of the impurity problem in ID systems have attracted renewd interest from the 
different branches of physics. |P, 0] The motivation of this paper is to investigate 

the impurity effect in the 1-dimensional interacting electrons under integrable boundary 
conditions. 

Generally, it is expected that at low temperature the quantum impurity Hamil- 
tonians are renormalized to critical points which correspond to conformal invariant 
boundary conditions. M, B] In fact, the impurity effects have been discussed from the 



analysis of the different boundary conditions. H, H, [TQ] During the last decade the 
Bethe ansatz techniques for the integrable open chains have been developed. ITTl , [T^ , 
131 , p^ , |15| In Ref. |^, the magnetization of an anisotropic Heisenberg model with 
open-boundary condition was derived. The result is generalized to the supersymmetric 
t-J model with open boundaries and the bulk and surface magnetizations are obtained. 



T0| We can also discuss the Bethe ansatz equations for the ID Hubbard model under 



some open-boundary conditions. [T^, |TB|, ^ However, the thermodynamic proper- 
ties of the t-J or Hubbard models have not yet been discussed under the open-boundary 
condition. It seems that under the open boundaries the free energy becomes divergent 
due to an infinite number of zero modes. 

Let us discuss one of the most different properties of the ID impurity effect. In 



the Kondo problem the Wilson ratio plays an important role. pO| Based on the local 
Fermi liquid theory, the local impurity effect in the 3 dimensional free electrons is fully 
characterized by the single parameter which is related to the ratio of the specific heat 
and the magnetic susceptibility due to the impurity. |^ For the ID Hubbard model, 
however, the local impurity effect could not be described by one parameter and it should 
be highly nontrivial. The low-energy spectrum is given by the Tomonaga-Luttinger 
liquid, which is completely different from the Fermi liquid. Furthermore, the impurity 
effect also depends on the Coulomb interaction among the conduction electrons. In 
this paper, we will evaluate the boundary effects to the magnetic susceptibility and 
specific heat of the ID Hubbard model under the open-boundary conditions, which 
should characterize the impurity effect in the ID interacting electronic system in the 
same way as the Wilson ratio for the Kondo problem. 
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Let us first review the Hamiltonian on the 1-dim. lattice with L sites [|TS 

L-l L L 

h ^ 

Here the symbols — /i and —"Pja (j = l,or L) correspond to the chemical potential and 
boundary fields, respectively. The symbol V denotes the Coulomb interaction. When 
= 'Pji (boundary chemical potential) or -p^^ = —pji (boundary magnetic field) for 
j = 1 and L, we can solve the above Hamiltonian by the Bethe ansatz method. The 
solution of electrons with M down-spins has wave numbers kj for j = 1, . . . N and 
rapidities Vm for m = 1 . . . M. The Bethe Ansatz equations for it are given in the 
following, ig 

(e*^jpil + l)(e"*''j +Plt) 



(sin kj -Vm + iU/4:){sm kj + Vm + iU/^) , . 

m= 



\^ (sin kj — Vm — iU / A) (sin kj + Vm — iU / A) 



(C+ - Vm - ^^/4)(C- -Vm- lU/^) fr (^m - Vn + lU /2){Vm + + ^^/2) 

SL {vm — sin kj + iU/A) {v^ + sin kj + iU/A) 
jj[ (f m — sin kj — ill/ A) {vm + sin kj — ill/ A) 

(3) 



where 



00 for pi| = pii (00 for = pLi 

for p,, = -p,, I - for p,, = -p,, 

(4) 



In Ref. [|T^], the Bethe ansatz equations (0) and (^ are systematically derived by using 



the refiection equations. 

2 The magnetic susceptiblity 

Now, we discuss the derivation of the magnetic susceptibility of the Hamiltonian (|l]) 
at zero temperature. We shall evaluate the boundary contributions to the suceptibility 



3 



using the method of Ref. where the susceptibihty is derived under the periodic 
boundary condition. For the evaluation of the susceptibihty we assume that the electron 
density n = N/L is less than half-filling (0 < < 1). 

Let us consider the ground states of the open-boundary Hubbard model for the 
repulsive {U > 0) and attractive (f/ < ) cases. The equations (0) and hold for 
positive and negative values of U. However, the ground-state solutions of the Bethe 
Ansatz equations are different for the two cases. For f/ > 0, we may consider only the 
case when the ground state is characterized by real k/s and real w^'s []• For U < 0, we 
may assume that the electrons may form singlet bound pairs in the ground state; the 
ground-state solution of the Bethe Ansatz equations consists of real k/s, real Vm^s and 
pairs of complex momenta k^ 



sm{k. 



(5) 



where u = |f/|/4. 

We solve the Bethe Ansatz equations based on the assumptions of the ground state. 
Hereafter, a subscript r = " > " ( r = " < " ) stands for the positive (negative) U 
case. For the positive U case, let p>,l(A;)i and p>,L(f )2 denote the densities of electron 
momenta kj and that of down-spin rapidities Vm , respectively. For the negative U 
case, we denote by p<,l(^)i and p<,L(f )2, the densities of real momenta kj and that of 
the string centers Vm, respectively. |^ We denote by p^^^ the vector of the densities 
Pr,L = {Pr,L{k)i, pr^L{v)2) foT T =>,<. We uow take the asymptotic expansion with 
respect to 1/L. Then, we have the following integral equation 

Pr,L{k,v) = plL{k,v)+J^,{kMk\v')p,^L{k\v'). (6) 

Here the initial values of the densities p^i{k,v) = p^^{k,v) + Tj./L are given by 



'>,oo 



{k,v) 



71 





r'i{k,v) 



(7) 



<,oo 



{k,v) 



-Re 

TT 



lU] 



J 



T'Ak,v) 



iQ<,oiv) 



(8) 



^For some values of the boundary fields, the Bethe Ansatz equations may have pure imaginary 
solutions (boundary bound states) for the ground state. 
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The definitions of d/dkPrfi{k) and d/ dvQrfi{v) will be given in eq. 
while for r =< they are given in the following 



2vriP<.o(A:) 



2(1 cos/c) 



2Pit(Pit + cos A;) 



+ 



(0) for r => g, 

2^ cos /c 



1 + Plt + cos A; 1 + + 2pi| cos /c sin A; + 

(t; — 2 



1 + Plt + 2pLT\/l - (t; - y^l - (f - m)^ 



Pil + y^l - (t; - m)^ 



(9) 



4m 



l+pf^ + 2pi|^l- 



f — 



f — 



+ 



+ 



2(m + <_ 



t;2 + 4u2 v'^ + {u + iC+f + + 

We note that the kernel was given in |]18| and = a^K^a^ [|2^. The parameters 
Q and B |]l8l for the upper- or lower-bounds of the integral intervals have the following 
constraints 



[Q 1 1 

/ P>,L{k)idk = 2n + — , / p>,L{v)2dv = n - 2s + — 

J-Q L J-B L 

/ P<,L{k)idk = As + — , I p<,L{v)2dv = n-2s + — 

J-Q L J-B L 



(10) 



where s has been defined by s = (A^ — 2M)/2L. The ground-state energy Er for 
r =>, < is given by 



% = ^[l-/^s-V2] + (e°,p>,J 
^ = j[l-ps-p + 2VlT^] + {e%p 



'111 



, t'T, 



where /i^ = /i/2 — /;,/4. The dressed energy satisfies 6^ = 6° + with the initial 
values 



jig — cos k 
h/2 



/is — COS 



-2i?e./l 



[V — lU] 



(12) 



By minimizing the ground state energy with respect to the variable s, we have the 
following functional relation between s and h through the constraints (|IOD. 



h 



det UQ,B) 



for 



r =>, < . 



(13) 



^d/dvQyfi{v) is given by d/dvQQ{v) of n = 1 in eq. (C 
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The dressed charge matrix is defined by the relation $,j.{k, v) = 1+Kj{k,v\k', v')$,^{k', v') 
and the symbols (rj are given by the matrix elements of the dressed charge as fol- 
lows; when r => (yj = (^>)ji for j = 1, 2, and when r =< C<j = (^<)ji + 2(,^<)j2 for 
j = 1,2. The symbols are defined by 




e° + 



h/4: 





(14) 



We can calculate the magnetization s through the equations ([TSD and (p^Ol) . We 
thus derive the magnetic susceptibility Xr,L of the finite-lattice of L sites, for the r => 
and r =< cases. Here L is a large but finite number. 



Xr,L 



dh dQ dh dB \ 
dQ'dS^dB'dSj 
2vr,i{Q) CAB) 



+ 



for 



PrAQ)i det^ B) Pr,L{B)2 det' ^,{Q, B) 



(15) 



The Fermi velocity Vrj are given by the j-th component of defined by 



.'0 



+ KrVr, where denote the vector whose first and second components are 
d/dker{k)i and d/dver{v)2, respectively. We note that if we specify the density n 
and magnetization s, then the parameters Q and B in eq. ([T5|) are defined by (10). 

Let us discuss the boundary contribution 6xr to the susceptibility for both the 
repuslive and attractive cases. We assume that the finite and the infinite systems have 
the same n and s. Then we may formally define Sxr by the following 



^Xr XrjL^QjB) Xr,oo(Qoo) -^oo)) 



for 



=>,< 



(16) 



Here Qoo and B^o are the interval parameters for the infinite system, which are given 
by eq. ( p!OD after taking the infinite limit: L ^ oo. For the case of nonzero magnetic 
field {Boo 7^ oo), we may evaluate 6B = B — B^o by taking the derivatives of eq. (p!0|). 
In terms of the dressed charge we have 



6Q 
6B 




( 



Tr{k)idk 



\ 



1 



J 



det^^ 



for 



1 - / Tr{v)2dv 
\ "'-Boo ' 

(17) 

Here the matrix elemnts of the dressed charge are evaluated at Qoo and B^o- For U > 0, 
by the Wiener-Hopf method we can show that under zero magnetic field B is as large 
as log L. We also note that for some values of the boundary fields, the magnetization 
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s of the finite system can take a nonzero value under zero magnetic field: h=0, in 
general. 

We find that the boundary contribution 6xr to the magnetic suceptibility contains 
both the charge and spin parts. We recall that the densities of the rapidities contain 
the 1/L-terms, which come from the open-boundary condidtion. We note that under 
the periodic boundary condition, the finite-size corrections of the densities do not have 
1/L-terms; the first nonzero order is given by 1/L^-terms. Therefore, the 1/L-term 
together with 6Q and 6B will reflect the effect of the open-boundary conditions. It 
seems that the result is different from the perturbative calculation of the 6x in Ref. 
0. However, the Sx in Ref. is obtained by using the bosonization method, where 
some limiting procedures are employed. Thus, it is not easy to point out the most 
important reason why they are different. We shall discuss this possible discrepancy in 
later publications. 



3 The specific heat 

In the rest of the paper, we show how to calculate the low-temperature specific heat of 
the open-boundary Hubbard model. We consider only the repulsive case {U > 0). For 
the negative U case, we can derive the similar results making use of the particle-hole 
transformation. Under the zero boundary fields case, the Bethe ansatz equations (2) 
and (3) are equivalent to those of the periodic case. Thus, the solutions of the Bethe 
ansatz equations have the same structure such as in the periodic Hubbard model. 



There are both real and complex solutions for the momentum kj and rapility Vm- They 
can be classified into three groups: real momenta k, n-A strings and n-A-k strings. The 
n-A string solution for rapility v is given by 

X^^' = A™ + m(n + 1 - 2j), j = 1, ■ ■ ■ , n. (18) 
The n-A-k string solutions for the momentum k and the rapility v are defined by 

■^m ^ = + + 1 - 2j), j = 1, ■ ■ ■ , n 

fc;;^2i+i ^ 7r-sin-^(A;, + m(n-2j)), 0<j<n-l 
klf^ = sm-\X'^ + tuin-2j)), I < J < n - 1, 
^m^" = TT — sin""^(A'^ — mn). 



(19) 



Here A„ and A^ are the centers of an n — A string and an n — X — k string, respectively. 
The symbols M„, and M' are the numbers of n-A strings, n-A-k strings and the all 
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A-k strings, respectively. The k/s form real momenta for j = 1, ■ ■ ■ , N — 2M' . We note 
M' = X^n^^n- recall u = U/A. In this paper, we only consider the case where all 
solutions are given by the three groups with the strings. 

Taking the asymptotic expansion ( up to the order of 1/L ), we can derive the 
following integral equations of the densities of particles and holes from the Bethe Ansatz 
equations (H) and @. 

p\k) = -p(A;) + - + -— Po(A;) + E / ^n(sinA;-A)cosA;K(A) + a„(A))rfA 

TT L dk n^-^ 

1 d r °° 

^n(^) = T7V^o(A) + / ^;(sinA; - X)p{k)dk - ^ A„„a^(A) 



^ \^l-(A-mn)2/ 

/oo 
6'^(sinA; - X)p{k)dk - ^ y4„„cr^(A) 
m=l 

PfM = ^ 1 (l+PiTe'^')(m + e^') 

oo 

- ^ e„(sin fc)(2 - 5(M^, 0) - 5(M„, 0)) 



m=l 

oo 



m,=l 



1 ^^^^ (A + C+ + ^uin - 2j))(A + C- + ^uin - 2j)) 



27l"^ (A - C+ - iu{n - 2j))(A + C- - - 2j)) 

oo 

Q^(A) = E0n-(A)(1-5«,O)) 

m=l 

1 2n + 



1 y^j^ (A + C+ + m(n - 2j))(A + C- + m(n - 2j)) 
2711 fr{ (A - C+ - m(n - 2j))(A + C- - tuin - 2j)) 
9 f III 

(20) 

where 6'„(x) = 2tan ^{x/nu)/27r and 6„m(a;) = (1 - Snra)0\n-m\{x) + 26\n-m\+2{x) + 
■ ■ ■ + 29n+m~2{x) + 9n+m{x). The symbol 6{j, k) denotes the Kronecker delta. We note 
that p, (Jn and 5"„ are the particle densities of the real momenta fcj's, the centers of the 
n — A strings, and the centers of the n — A — k strings, respectively; p^, and 5"^ are 
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the hole densities of them, respectively. It is remarked that in the derivation of the 
Bethe ansatz equations have assumed that the string solutions for the finite 

system could have small deviations from those of the inifnite system given in (IS). 
The total energy E of the system is given by 

- = -i^{l-f,,){l-6{N,2M'))--j:n{l-6{M^,0))j 



Y.i2Jl + (nuy + n^^){l-5iM:^,0)) 



n=l 



{fis — COS k)p{k)dk + — "^n / (T„(A)(iA 

-TT 2 J -00 

\Jl — {\ — inuy + n/j^ an{\)d\ 



(21) 



+ E 

n=l 



It should be emphasized that the sums for the zero modes in ([21|) do not become infinite 
because we take L a large but finite number. 

Minimizing the thermodynamic potential Q = E — TS, we can show the following 
thermal Bethe Ansatz equations for the ratios of the particle and hole densities ( = 

P^/p, r]n = (Tn/cTn, and fin = O-n/^n 



lnC(A;) 
ln?7i(A) 

Inf^i(A) 
ln?7„(A) 
In r7n( A) 

lnr7„(A) 

In fyn( A) 



2cosfc 7r(A-sinfc) / l+fji 4:Re{^ 1 - {X - iuy \ dX 



T 



+ 



sech( 



2u 



In 



T 



Au 



s* i^n{l + r]2{X)) - J \n{l + C\k) cosk6{X- sin k)dk 

s* (ln{l + fi2iX)) - ln{l + C (k) cos k 5 {X~ sin k)dk^ 
s * (ln(l + Vn+iiX)) + ln(l + Cn^i(A))) 
. * (ln(l + Vn+i{X)) + ln(l + C„_i(A)) 



h 
n— 
T 



n- 



4:U + 2p 



(22) 

where the convolution s* is given by s* f{x) = J {sech.{7i{x — x') /2u) / 4u) f {x')dx' . The 
last two limits in (E^) should be understood as the asymptotic limit of n consistent 
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that of the notation n — oo means that n should be taken a very large but finite 
number. We should note that the variation of the particle and hole densities should be 
taken over all the positive values of rapidities since they are even functions; we have 
such as p{k) = p{—k), and so on: a{v) = cr(— 1>), and a{v) = a{—v), p^{k) = p^{—k), 
a^(v) = a''{-v), and a''{v) = a''{-v). 

Substituting (pO|) into Q we have the asymptotic expansion of the thermodynamic 
potential ul = Q/L with respect to 1/L 

1 / h °° \ 

= -(^(l-/x,)(l-5(iV,2M'))--E(l-'5(M„,0))j 

1 °° I 

-- 5:(2v'l + (n«)2 + n/i)(l - 5(M;, 0)) 

^ n=l 

(23) 

We now introduce k = Tln^ and ei = Tlnr^i. We denote the zero-temperature 
limits of K and ei by k^^^ and e^°\ respectively. Hereafter we assume h ^ T. We can 
calculate the specific heat Cl through Cl = —Td'^UL/dT'^. We give the final results: 
a: -M < -2 - h/2 



ujl{T, h,p) = ujL{0,h,fi) 
where 
■5/ = 



a 4) f in (l + e-" exp{i±i!^}) (24) 



TT 



1/(1+ Pit) , Plt ^ -1 1 _ J Pit/(1 +Pit) , Pit ^ -1 
1/2 , m = -i J 1 1/2 , PiT = -l 



Let Coo and (5C denote the bulk and the boundary specific-heats, respectively. Then 
from (|2^) , we see the following 

Cl = C^ + SC = C^l + ^Sf). (25) 



b: e5°\0) > 0, -^i>-2-h/2 
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In this region, the thermodynamic potential is given by 



LJL{T,h,fi) = ujL{0,h,fi)- 



1 d 



12 sinQ \ TT L dk 



PoiQ) 



g{0) + -6g{0) 
2u 



In 1 + e""" e~ 



e(°)(0)/TA 



dx 



dk 



2tx J-q (A - sinfc)2 + ^2 



2u 



2J-Q{X-smky + u^dk ^''2n 2dX 



(26) 

Here the paramter Q is the zero of K^^\k). From this expression, we find the ratio 

5C/Coo 



C(X) dk 

c: h > 4(V1 + - u), -n>2-h/2 

In this region, we find 



(27) 



UL{T,h,fi) 



u;l{0, h, /i) - T'/'n-\l + Lihl^] 



rin{l + e"e-''')dx 
Jo 



-T3/24(1 + m2)1/V-1(i + £) r ln(l + e^e-''')dx 

L Jo 



a 



2 + u-h/2 



T 



TT 



h - A{y/TTv^ - u) 
~ f 

2« iPoik),, , d^. 



(2J 



TT sin k + 2n 



'^^'^^dk+^Qli^) 



dX 



d: e'i\0) < 0, KW(7r) > 

Let us denote by B the zero of ei'\X). Then from the thermal Bethe-Ansatz 
equations (^3|) , we obtain 



u;l(T, h,fi) = ljl{,0, h,fj,) 



pUQ) 



(29) 



where the density functions are given by p1{k) = p>,l(/i;)i and p1{v) = p>,L(f )2, where 
Py = {p>,L{k)i, p^^l{v)2) is the density for the ground state of the repulsive case. 
(See also Ref. ||). ' 
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e: 4V1 + -u> h:^T, k^^^tt) < 

In this region, the free energy can be evaluated as 



(30) 



where po and af^^ are defined by 



a{°)(A) = ao(X)+ f R(X - X')a[''\x')dX 

J\X'\>B 

^e(o)(A) = r sisin k-X) cos kdk+ f RiX - X')^/^\x')dX' 

J —IT J\x'\>B dX 

<Jo(A) 



U{X - sin k)dk + i ( J^Qi(A) + it: * J^Qo(A)) 



Po(/^) 



+ 



1 / d 



TT L \dk 



Po{k) - Qo{k) 



i-oo / I 

+ COS k / ai(sin A; — A) I (Jo (A) — - 

J —OO \ Li 



dX 
(31) 

and Q{k) = cos A; / s(sinA; — x)d/ dxQ\{x)dx. Here / * g{x) — J f{x — x')g{x')dx' , 
s{x) = scch{nx/2u) / Au and R{x) = s * 2u/ (27r(x^ + u"^))- 

From the low-temperature expansion of the thermodynamic potential, we find the 
boundary contribution to the specific heat at low temperature 



5C 



1 

L 



case a 
case b 

case c 



(32) 



+ 



/{C + C) , cased 



case e 
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Here Coo denotes the bulk specific- heat, and and denote the contributions to 
the bulk specific- heat Coo from the charge and spin parts, respectively; Coo = C'^ + C*. 
The symbols (Tj°]o(i?) and 6ai'\B) are the bulk and the 1/L parts of ai^\B). 

The boundary contributions to the magnetic susceptibility and the specific heat 
depend on the boundary fields and the electron density. In regions a, and c, the 
specific heat is proportional to T~^/^e~"/"^ where a is a positive constant. In regions 
b, d and e, the specific heat depends linearly on temperature. For the boundaries of 
the regions between a and 6, h and c, h and d, c and e, and d and e, the specific heat is 
proportional to T^/^. We recall that for some regions of the boundary fields, there may 
exist other types of solutions (boundary string states) of the Bethe ansatz equations. 
P^ , ^ We can calculate the contribution from the boundary string states 

simply by modifying the term t° in our derivation. 

In summary, we have studied the boundary contributions to the magnetic suscept- 
bility and the specific heat for the 1-dimensional Hubbard model under the general 
open-boundary conditions. They are calculated analytically and explicitly. From the 
results we can discuss exactly the impurity effect in the 1-dimensional Hubbard model 
or in the interacting electrons in ID. 

We are grateful to Dr. K. Kusakaba for helpful discussion on the boundary bound 
states. R. Yue was granted by the JSPS foundation and the Monbusho Grant-in-Aid 
of Japanese Government. 
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